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Abstract. The purpose of this article is to explore the optimal choice of shape parameter which 
is an important and longstanding problem in the theory of radial basis functions(RBF). We already 
handled it for multiquadric and Gaussian in [9j [10l [11] [T2l US] . Here we focus on shifted surface 
spline and present concrete criteria for the choice of shape parameter. 
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1 Introduction 

The so-called shifted surface spline is defined by 

h(x) := (-l) m (|a:| 2 + c 2 )^log{\x\ 2 + c 2 )^, A e Z+, m = 1 + ^, c> 0, x € R n , A, n even (1) 
, and 

h(x) := (-1)^-^ {\x\ 2 + c 2 ) x -% , Xe Z+ = {1,2,3, • • •}, n odd and A > - (2) 

, where |x| is the Euclidean norm of x, and A, c are constants. The number [A — denotes the 
smallest integer greater than or equal to A — 

This function was introduced by Dyn, Levin and Rippa in [3] for x £ R 2 . Then it's extended to 
R n for n > 1. For further details we refer the reader to [2 El [5j [T8J [19j [111 [22l [23l [24] . 

Note that (2) is just multiquadric and will not be explored here. In this paper we only deal with 
(1). As is well known, h(x) is conditionally positive definite of order m = 1 + 4 where A is of our 
choice. The constant c is just the so-called shape parameter whose optimal choice is of our primary 
concern. 

Since h(x) is c. p. d. (conditionally positive definite) of order m, there exists a unique function s(x) 
of the form 

N 

s(x) — p(x) + Cjh(x — Xj) (3) 
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which interpolates any scattered set of data points (xi, y±), • • • , (xn^n) as long as {x\, • • • , xn} is 
polynomial-nondegenerate in R n . Here p(x) € P m — l is a polynomial of degree < m — 1 and c^-s are 
chosen so that 

N 

^2cjq(xj) = 
i=i 

for all polynomials q in P m _i. 

We will show that whenever x and the data points are fixed, the upper bound of \f{x) — s(x)\ 
will be made minimal by a good choice of c. Here f(x) denotes the underlying function of the data 
points. 



1.1 Basic Ingredients 

Definition 1.1 For n = 1, 2, • • • , the integers 7„ are defined by 71 = 2 and 7„ = 2n(l + 7 n _i) /or 
n > 1. 

The numbers 7„ will appear in our criteria. 

In this paper all approximated functions belong to a semi-Hilbert space Ch. m (R n ), abbreviated 
as Ch,m, where m denotes the order of conditional positive definiteness of the function h. There is a 
semi-norm on Ch, m , denoted by || • \\h- The definition and characterization of this space can be found 
in [HI [16]. We will not repeat them. A subspace of Ch, m , denoted by B a , is of our special 
concern and is defined as follows. 

Definition 1.2 For any a > 0, 

B a := {/ e L 2 (R n ) : /(£) - for all |£| > a} 

is the set of band-limited functions, where f denotes the Fourier transform of f . 

We need some constants as follows. 

Definition 1.3 Let n,X and m be as in (1). The constants p and Ao are defined according to the 
following three cases. 

(a) Suppose n — A > 3. Let s — |~ "~^~ 3 ] . Then 

s (2m + 2 + s)(2m + l + s)---(2m + 3) 
P=1+ 2rJT3 and A ° = ■ 



(b) Suppose n - A < 1. Let s = - 1 1 ^^-] . Th 



en 



p = 1 and Ao = 



(2m + 2)(2m + 1) • • ■ (2m - s + 3) ' 
(c) Suppose 1 < n — A < 3. Then 

p = 1 and Ao = 1. 

Finally, the Fourier transform [6] of (1) should also be introduced. It's of the form 

m = i{\nm- x - n K, n± ±{c\a\) (4) 

2 

where l(X, n) = (2n)-% ■ 2? ■ (f !) 16 and JC v (t) = t v JC v (t), K v {t) being the modified Bessel function 
of the second kind pQ . 
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2 Fundamental Theory 



The cornerstone of our theory is the exponential- type error bound raised by the author in |14j . We 
cite it directly with a slight modification to make it easier to understand. 
In the following theorem 

d — d{E, X) := sup inf \y — x\ 

denotes fill distance. 

Theorem 2.1 Suppose h is defined as in (1). Then, given any positive number bo, there are positive 
constants dg and lj, < u> < 1, which depend on b , for which the following is: If f € Ch^ m and s 
is the h spline defined in (3) that interpolates f on a subset X of R n , then 

\f(x)-s(x)\ < y/l(X,n) • (2tt)* • V^T- V^' ^ ■ lj* ■ \\f\\ h (5) 

holds for all x in a cube E C il, where £1 is a set which can be expressed as the union of rotations 
and translations of a fixed cube of side b^, provided that (a)E has side b > bo and (b)0 < d < d<j. 
Here, a n denotes the volume of the unit ball in R n , c was defined in (1) and Ao was defined in 
Definition 1.3. 

The numbers do and ui can be expressed specifically as 

d ° = 6C ln (m + iy 

where 

C = max|2p'^e 2 ^", J-}, p' = L 

The numbers j n and p were defined in Definitionl.l and 1.3, respectively, and m = 1 + -| as in (1). 

Remark: (a) The seemingly complicated theorem is in fact not difficult to understand. We expressed 
each constant specifically for later use. (b)The set X of centers is usually contained in the cube 
E. However it's harmless to extend it to a more general form as a subset of R n . (c)\\f\\h is the 
semi- norm of / mentioned in subsectionl.l. 

On the right-hand side of the inequality (5), the only things dependent of the shape parameter 
c are ca, uii and \\f\\h- Also, do is dependent of c. It's not hard to imagine a good choice of c 
will minimize the error bound (5). However it's nontrivial to express this error bound by an explicit 
function of c, and clarify the relation between do and c. Further treatment of (5) is necessary. 

Lemma 2.2 Let a > and h be as in (1). For any f £ B a , 

i 

ll/IU < C (m,n) ' (^) 4 ' ' ■ • (6) 

, where 

Co (m, n) 




(27r)~"Vm! 
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Proof. By Corollary3.3 of [TB], 



1/2 



|a| — m 



i 2 (p) 



(-)-{ E 5/ 



1/2 



1/2 



C 2Q l/(OI : 



< 



(27r)-"Vm! 



|C| 2 ™/(A,n)|C|- A -"/C^( C |f|) 



(27r)-™Vm! 



1/2 



iei- A -"/c^( c |^i) 

l/(OI 2 



< Co(m,n) C ^|y |/(£)| 2 |£l 
(27r)-»V^i 

C {m,n) - 



iei- A -"Wci)^/c^(ciei) 
i 



1/2 



1/2 



2 ' v^l«l 



■c?0 by [20] where 



Cq (to, n) 



1/2 



If / € -B CT , then 



|/IU < C (m,n) 



a 4 c * e ' ||/|| L 2. 



Corollary 2.3 Lei cr > 0. If f £ B a , then (5) can be transformed into 

\f{x) - s{x)\ < v^(27r)~"V2na„cr i± T ±A \/A^c i ^ ±A e^w3||/|| L 2. (7) 

Proof. B a C Ch m is a simple result of Corollary3.3 of [16] . Now (7) is just a combination of (5) and 
(6). tt 

On the right-hand side of (7), the only things dependent of c are c ^~ , e~ and w^, where d 
denotes fill distance. It seems that one can abstract a function of c from the long expression of (7) 
after every thing independent of c is fixed, including d. 

However as c changes, the upper bound do of d changes also. As required by Theorem2.1, d < do 
where 

d ° = 6C 7n (m + l)- 
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The number C — » oo, and hence d —¥ 0, as c — > + . In order to keep d < d , the minimal possible 
choice of c is 

c := 12p^/ne 2 " 7 "7„(m + l)d (8) 

where d is fixed. Therefore, we can choose c from the interval [co, oo) only. 

There is a logical problem about d, do, c and bo- Before c and 60 are given, do is unknown and we 
do not know whether d is eligible. This question is not difficult to resolve. For any bo > 0, we require 
d < 6C 7 „(m+i) wn ere we let C = ^ temporarily. Then d < iry J^ +1 y Note that 2p'^/ne 2n ^ n = ^ 
iff c = ZboP\fne 2n ' 1n . With the same b and c > ZboP\fne 2nin , we have 



C* = max|2p'V^e 2 ™ 7 " 

and 

rf 



2 12 



4 7 „(m + l) 

Then d < d is satisfied. 

With the same b and c < 3bopy/ne 2njn , we have 



C = max ^2p'^e 2n ^, A J = Ipfyfa 



ne 2 ™ 7 " 



and 

d 



6C7„(m + I) I2p^ie 2n ^ ln (m + iy 
In order to keep d < do, we require 

c > c := I2/V^e 2 " 7 "7„(TO + l)d 

where d was given above satisfying d < 47 (m+i) ■ Therefore, once the cube side bo is given, we first 
choose d < 4j ^+1) ■ Then put the restriction c G [co, 00) where Co = 12pv / ne 2Tl7 "7„(m + l)d. The 
condition d < do will always be satisfied. 

Theoretically d, and hence co, can be arbitrarily small. However practically the problem of ill- 
conditioning when constructing s(x) has to be considered. In this paper we explore theoretically the 
optimal choice of c for c < c < 00. 



3 Criteria of Choosing c 

In the long expression of (7) the only things dependent of c are c « +A , and . Let's define 

MiVfcj^c^.e^'W' (9) 

and call it an MN function. The optimal choice of c will then be the number minimizing MN(c), 
when every thing independent of c is fixed. 

The value of u highly depends on the cube side bo- Theoretically, uj — > and is very influential 
as c — > 00 and 60 - > 00 simultaneously. However, a lot of time bo is fixed and cannot approach 00. 
Therefore we divide our criteria into two classes. 
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3.1 b fixed 

Recall that C = max{2p' \/ne 2nin , -^} and p' = |. The two values 2//^ne 2 " 7 ™ and ^ are equal 
if and only if 

c = S&o/V^e 2 " 7 " =: ci (10) 
Then, for c e [cq,Ci), C = 2p' y /ne 2n ~ fn . For c <G [ci,oo), C = g|-. So we have 



3fco ' 

(|) 12pV5Te^T"-y„d if C g[c ,Cl), 

. ifce[ci,oo) 
e ^^"in^j ifce[co,ci), 
e if c e [ci, oo). 



Putting this result into MN(c), we thus have 



MN(c) = { c 4 ' e V 2 7 if ce [co.ci), (n) 

c 4 . e 2 • e 4 ^" d irce[ci,ooj 



which is a continuous function. Our goal is to find c which minimizes MN(c). 

As for the fill distance d, we require d < - A — , b ° ... once the cube side b is given, as mentioned 
in the end of section2. 

We have the following cases where n, A, m, cr, 7„, p and (i were defined in (1), Definition!.. 1, 
1.2, 1.3 and Theorem 2.1. Moreover we introduce a constant 

cr Zn| 



2 12pV^e 2 "T-7„d' 



Casel. 



n + A > and fc > Let / e _B CT be the approximated function, b > be the 

be the fill distance. If 1 — n + A>0 and fc > 0, the 



cube side as in Theorem2.1, and d < 



b 



47„(ro+l) 



optimal choice of c for c e [c , oo) is c = 12py/ne 2nirl j n (m + l)d 



Reason: In this case MN(c) is an increasing function in its domain [co,oo). We therefore choose 
c = co, defined in (8), which minimizes MN(c) in (11) and the error bound in (7). 



Case2. 



n + A > and k < Let / e B a be the approximated function, b > be the 

+ A > and k < 0, the 



cube side as in Theorem2.1, and d < -. — fa . . , be the fill distance. If 1 — 
optimal choice of c for c e [co,oo) is the value minimizing MN(c) over the interval [co,ci] where a 
was defined in (10). 

Reason. In this case MN(c) is increasing on [ci,oo). Therefore the minimum value of MN(c) 
in [cq,oo) must happen in [c ,ci]. 
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Remark: (a)This case rarely happens because k is usually positive. (b)In fact, the optimal c in 
Case2 can be obtained exactly. Let g(c) := c ^~e ck which is just MN(c) on [co,ci]. Then 5(c) 
is increasing on (0, 1 ~^ A ] and decreasing on [ A , oo). So the optimal c is (i)c if c\ < 1 ~ T ^ A , 
(ii)ci if < c , and (iii)ci if g(a) < g{c ), and c if g(co) < g{c\), for c < < d- 



Case3. 



1 - n + A < 



and k = Let f € B a be the approximated function, b > be the 
cube side as in Theorem2.1, and d < 4j (^+1) be the fill distance. If 1 — n + A<0 and k = 0, the 
optimal choice of c for c e [c , oo) is the value minimizing c 4 +A e ^f" on the interval [ci, oo). 

Reason: In this case MN{c) is decreasing on [co,ci). Therefore the minimum value of MN(c) 
happens on the interval [ci,oo) where the essential part of MN(c) is c 4 +A e ^f\ 

Remark: (a) Although k rarely equals zero, we can make it zero by choosing a in an appropriate 
way. It will make things easier for the optimal choice of c. (b)If 1 — n + A < 0, it can be shown that 



9(c) :-- 



~e 2 is decreasing on (0, 



-l+n-A] 



and increasing on 



-l+n-A 



2cr j^^v^t,^! 2(T 

Therefore the optimal choice of c for c e [co, oo) in Case3 is in fact 



oo), with g' {■ 



-l+n-A 
2a 



)=0. 



max 



-l+n-A 

2cr 



,Cl} 



Case4. 



1 - n + A < 



and 



k > 



Let / e B„ be the approximated function, 6 > be the 
jhrn t" 3 tnc nu distance. Let 51(c) := M7V(c)| ce [ COiCl ) and 



cube side as in Theorcm2.1, and d < 7— r 

52(c) := MiV(c)| c£ [ Cl;(X) ). Ifl-n + A<0 and k > 0, then 

« SK 1 ^) = 0, 5i(c) < for c e (0, ±=^), and ^(c) > for c e (^^,00), and 

(ii) .92(^7^) = 0. 92(c) < for c e (0, and .g 2 (c) > for c e 00). 

In this case the minimum value of 51(c) on [co, Ci) happens at (a)c = co if 1 ^""^ A < Co, (b)c = 1 ~"^ A 
if c < 1 ~"fc A < c i ; and (c)c = Ci if Ci < 1 ~^ A . Also, the minimum value of 52(c) on [ci,oo) 



-l+n-A 



< ci, and (b)c = if ci < 



happens at (a)c = ci if — ^ . u i7 — — ^ 

Let c* e [co,Ci) minimize 51(c) and c** € [ci,oo) minimize 52(c). Then the optimal choice of 
c e [c ,oo) is (a)c* if 51 (c*) < 52(0**), and (b)c** if 5 2 (c**) < 51 (c*). 

Reason: By direct differentiation, we get 



g' 1 (c)=c * 



n + A 



+ ck 



and 



5 2 (c) = c" 



_1 CCT _^ 3 

e 2 e 4 ^« d 



1-n + X 



+ 



ca 



4 2 

The seemingly complicated criterion is then just a simple result of the two derivatives. 



Case5. 



1 -n + A < 



and 



k < 



Let / € -B<r be the approximated function, 60 > be the 



cube side as in Theorem2.1, and d < 



bo 

4 7 „(m+l) 



be the fill distance. If 1 — 



optimal choice of c for c e [cq, 00) is (a)c = C! if 



-l + Tt-A 

2a 



< ci, and (b)c = ~ 1 +"~ A if d < 



A < and fc < 0, the 



-l+n-A 
2ct 



Reason: In this case MN(c) is decreasing on [co, ci) and its minimum value happens at a number in 
[ci, 00) which minimizes 52 of Case4. Our criterion thus follows immediately from Case4. jj 
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3.2 b not fixed 



In Theorem2.1 one can easily find that if 60 is not fixed and can be made arbitrarily large, C will 
become arbitrarily small by letting c — > 00. This will make ui — > and do — > 00, a very beneficial 
situation. In fact some domain Q does allow bo —¥ 00. For example fl = R n or 

O = {(xi, • • • , x n )\ Xi > for i = 1,2, • • • ,n}. 

This kind of domain is called dilation-invariant by Madych in [17] . For this kind of domain we 
have a different set of criteria of choosing c. 

By increasing the cube side bo, one can always keep C = 2// v / ne 2 ™ 7 ™ in Theorem2.1. Thus 



MN{c) 



l-n+A C 

c 1 e 



for all c£ [co, 00). 

Note that we never decrease bo because it will make the error bound (7) worse. 

As for the choice of d, there is no restriction. Once d > is given, Co will be fixed. For any 
c e [co,oo), the condition d < do in Theorem2.1 will be satisfied. However the smaller d is, the 
larger [cq, 00) is, making the criteria more meaningful. 

Once c is chosen, we let bo = 3pv ^,2n 7 „ such that 2 j o'y / ne 2 " 7 ™ = in the definition of C . 
Further increasing bo is not to be expected because more data points will be involved to keep the 
fill distance d fixed. 

Note that by letting k := | + 12p ^ e |„ 7rid , we have a more simple expression for MN(c) as 

MN(c) = c^^e ck . 
The way of choosing c is then divided in the following cases. 



Casel. 



1 -n + A > 



and 



k > 



Let / € B a be the approximated function, and the domain 
in Theorem2.1 be dilation- invariant. For any fill distance d > 0, ifl — n + A>0 and k > 0, the 
optimal choice of c for c € [co, 00) is c = Co. 

Reason: In this case MN(c) is increasing on [cq, 00). jj 



Case2. 



1 - n + A > 



and 



fe < Let / G Bp be the approximated function, and the domain 



fl in Theorem2.1 be dilation- invariant. For any fill distance d > 0, ifl — n + A>0 and k < 0, the 
larger c is, the better it is. 



Reason: In this situation, 



MN'(c) 



1 L e ck 



1 -n + A 



iff 



1 — n + A 
^4fc ' 



The function MN(c) is increasing and decreasing, respectively, depending on c < 

n± 
-4fc 



1-n+A 
4fe 



or c > 



1 n +\ Note that MN(cq) is a finite value and MN(c) — > as c — ► 00. Our criterion thus follows. 
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Remark: (a)Usually k > 0. Hence Case2 rarely happens. (b)If cq < 
small enough to be accepted. If it happens, c = c is also a good choice. 



l-n+X 
-4k ' 



MN(c ) may be 



Case3. 



1 -n + A < 



and 



k > 



Let / G B a be the approximated function, and the domain 
il in Thcorcm2.1 be dilation- invariant. For any fill distance d > 0, ifl — ra + A < and fe > 0, the 
optimal choice of c for c e [co, oo) is (a)c = x ~_ n ^ x if c < 1 ~"^ A and (b)c = c if c > 1 



-rt+A 
-4fc 



Reason: By the structure of MN'(c), one sees easily MN{c) is decreasing on (0, 1 _"fc A ] and in- 



creasing on _4fc A ; °°)- O ur criterion thus follows immediately. 



Case4: 



1 - n + A < 



and 



k < 



Let / € be the approximated function, and the domain 
il in Thcorem2.1 be dilation- invariant. For any fill distance d > 0, ifl — n + A < and k < 0, the 
larger c is, the better it is. 



Reason: In this case MN(c) is decreasing on (0, oo). Moreover, MN(c) ->0asc^oo. jj 

Remark: (a)Casel-4 are all based on Theorem2.1. Although c and bo appear before d in Theo- 
rem2.1, in the four cases d is given before c and bo are chosen. Once d is fixed, we chose c according 
to the four criteria. After c is chosen, we let bo = 3p ^i e 2 n -, n -(b)In order to keep the fill distance 
d fixed, more data points have to be added when the cube side bo increases. (c)Although in this 
paper we only explore c over the interval [co, oo), the curve of MN(c) can be used as a conjecture 
for c e (0, c ]. It probably holds even if we cannot prove it. 
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